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Abstract. In this letter we present the soft X-ray (5-35 A) spec- 
trum of the supernova remnant (SNR) IE 0102-72.3 in the 
Small Magellanic Cloud, acquired by the reflection grating spec- 
trometer (RGS) aboard ESA’s XMM-Newton Observatory. This 
extended-source X-ray spectrum of unprecedented spectral res- 
olution (A/ A A ~ 300) permits, for the first time, unabigu- 
ous identification and measurement of isolated emission lines 
and line complexes alike. The diagnostic power of performing 
spectroscopy using groups of emission lines from single ions is 
exemplified. In particular, the bright Lyman and helium series 
lines for light elements (Cvi, Ovil, Ovm, Neix, Nex and 
possibly Mg Xl & Mg XII) show peculiar ratios, where the val- 
ues [Is — np] / [Is — (n + l)p] are systematically weaker than 
expected for electron impact excitation. These measured ratios 
resemble signatures of recombining or charge exchanging plas- 
mas. We argue that charge exchange, given its large cross sec- 
tion and evidence for inhomogenious media within the SNR, is 
a likely mechanism for the observed emission. Also, the well 
known temperature diagnostics G(T C ) = (i + /)/r ofhelium- 
like triplets (Ovil & Neix) indicate high temperatures, well 
above the maximum emission temperature T m for each ion, 
and consistent with a purely ionizing plasma. The density di- 
agnostics R(n c ) = //i meanwhile, are consistent with the low 
density limit, as expected. 
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1. Introduction 

The supernova remnant IE 0102-72.3 (E0102) is a bright, oxygen- 
rich remnant located in the Small Magellanic Cloud (SMC), ap- 
proximately 60 kpc distant. Although it is young (~ 1000 years 
old) and clearly very energetic, it has no bright emission lines 
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of Fe K nor significant continuum above 5 keV. It has received 
substantial attention recently, with the availability of high qual- 
ity Chandra data (Gaetz et al., 2000; Hughes et al., 2000) as 
well as Hubble Space Telescope data (Blair et al., 2000). Its 
ionization structure in terms of spatial stratification was clearly 
seen for the first time with Chandra (Gaetz et al., 2000), and 
details of the outgoing blast wave’s kinematics and emission 
spectrum were used to argue that a large fraction of E0102*s 
energy is spent toward cosmic ray acceleration (Hughes et al., 
2000). It was also observed with the High Energy Transmission 
Grating Spectrometer aboard Chandra to explore the spatially 
resolved, line of sight velocity structures to isolated emission 
lines in the spectrum, which help to elucidate the initial ge- 
ometry and mass distribution of the initial supernova event, as 
well as probe the kinetics of various spectrally resolved compo- 
nents making up the SNR structure (Flanigan, 2000; Canizares, 
2000 ). 

In our work, we concentrate on the high resolution X-ray 
spectroscopy afforded by the RGS (den Herder et al., 2000) 
aboard XMM-Newton (Jansen et al., 2000), which performs with 
nearly nominal spectral resolution, even for extended objects 
such as E0102: We observe using the RGS to address the emis- 
sion mechanism of the SNR. 


2. Observation and Data Analysis 

E0102 was observed early in XMM-Newton '$ Calibration phase, 
during revolution 0065 (16 April, 2000) for a total of 42.8 ksec. 
All three EPIC (Turner et al., 2000; Struder et al., 2000) instru- 
ments and both RGS (den Herder et al., 2000) instruments were 
operated simultaneously. Here, our primary focus is on the high 
resolution X-ray spectrum of E0102, and we restrict attention 
to the RGS data 1 . 

i- The data were processed using custom software, originally 
developed for the analysis of RGS ground calibration data, which 
is nearly identical in function to the RGS branch of the Science 
Analysis System (SAS). Telemetered CCD events were read in, 

1 A separate letter on the EPIC observation of E0102 is included in 
this volume (Sasaki et al., 2000). 
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frame by frame for each CCD node, and were offset corrected 
on a pixel by pixel basis using median readout maps, compiled 
from about 40 DIAGNOSTIC (den Herder et al., 2000) images 
per CCD chip. This process nearly eliminates flickering pixels 
from the dataset. Gain and CTI corrections were performed to 
align the signal/energy scale across all CCD readouts. Then, 
event reconstruction was performed on connected pixels con- 
taining significant signal, and the composite event signals were 
calculated by summing up signals from individual pixels. 

The standard event grade combinations used were those 
which fall within a 2 x 2 pixel region, where two pixels diago- 
nally opposed to one another within the 2 x 2 were not grouped 
into composites. The event coordinates were then mapped into 
focal plane angular coordinates (dispersion and cross-dispersion) 
Using the star tracker attitude history updates for the revolu- 
tion, the coordinates of the source, and a prelminary boresight 
axis, aspect corrections were applied to the focalplane event 
angular coordinates. 3 

Events were extracted in the dispersion-pulseheight plane 
with masks utilized by the response matrix generator. The sur- 
viving events were then windowed in the focal plane, where 



channel wavelength (Angstrom) 


Fig. 1. The first order (m=-l) RGS spectrum of E0102. The data ob- 
tained from the two RGS instruments are plotted separately (RGS 1 in 
black, RGS2 in red). The 20-24A gap in the RGS2 spectrum is due 
to failed electronics for CCD 4 in that instrument. For clarity, the data 
are provided in both linear and logarithmic scales. Identifications of 
the principle lines are provided. 


background subtraction was performed using background sam- The w j n d 0 wed background contribution to the (m=- 1) total 
pling regions off to one side of the illumination pattern in the countrate for this extraction was nearly fiat, at 2.3 x 10~ 5 ct s” 1 chan” 1 
cross-dispersion direction. 4 The size of the spatial extraction OJ ^ equivalently, 2.3 x 10" 3 ct s” 1 A -1 at 15A 5 . 
was 58" in the cross-dispersion direction, which comfortably 


includes the entire spectrum. A spectrum file was created by 


computing the background corrected countrates in each disper- 3 The Soft x-Ray Spectrum of IE 0102-72.3 
sion channel. Finally, a response generator was used to pro- 


duce an observation specific response matrix that provided an 
array of nominal wavelength/energy values corresponding to 
each channel in the spectrum file. Line fluxes were estimated 
by comparing an emission model, folded through the response 
matrix, directly to the spectrum in the usual manner. This pro- 
cedure was performed for first (m=- 1 ) and second (m--2) or- 
ders for each spectrometer. 

To reduce the contamination from sporadic background rate 
fluctuations in which the RGS counting rate doubles from the 
quiescent rate of about 1 ct s _1 RGS -1 , a set of good time in- 
tervals (GTIs) were chosen to filter out the high background 
data. This reduced the effective exposure time from 37.9 to 29.7 
ksec for both RGS instruments. 


Using the analysis procedure described above, we generated 
the RGS spectra shown in Figures 1 and 2. The spectrum is 
dominated by a small number of emission lines, the brightest 
of which are resonance transitions of highly stripped light ions. 
Nitrogen, which is cosmically less abundant than oxygen by 
a factor of a few, is is absent in our spectrum (N VII Lya is 
expected at 24.7 A). Only very weak Fe L lines are detected, 
which are dwarved by the strong transitions of the light, even 
Z elements. We detect bright lines from the Lyman series of 
C vi, O vill, NeX and MgXll, and from the helium series of 
O vn, Ne IX, Mg XI and Si XIII. Most of these ions are repre- 
sented by several emission lines in the spectrum, clearly re- 
solved and measurable in this phenomenally clean spectrum. 
Single ion spectroscopy permits independent measures of as- 
trophysical quantities that affect line production, such as elec- 
tron temperature and density, or radiative transport within the 


2 The dispersion coordinates are based on the known geometry of 
the RGS within XMM-Newton (den Herder et al., 2000). 

3 Based on attitude history updates of the star tracker system, the 
spacecraft pointing was extremely steady, drifting only 0.1" and 0.2" 
(rms) in the dispersion and cross-dispersion axes, respectively, over 
the observation. The respective peak-to-peak attitude swings were 
only 0.5" and 0.7". 

4 The target was ~ 1' off-axis in cross-dispersion, so the spectrum 
illuminated the RGS focalplane camera on one side of the CCD array. 
Consequently, background sampling was performed in an asymmetric 

manner. 


object and in any intervening material. 

To synthesize the appropriate RGS response for E0102, we 
need to quantify the angular distribution of the X-ray emission, 
which is convolved with the other quantities contributing to the 
line spread function. As a preliminary approximation, we used 
a public Chandra image (Obslds 1231 & 1423) as an energy in- 
dependent distribution, characterized by ~20" HEW and ~33" 

5 Also equivalent to 1.4 x 10" 4 cts" 1 per resolution element. This 
is higher than the typical background level because of the enlarged 
extraction volume for this target. 
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RGS channel wavelength [A] 

Fig. 2. Detail of the 8-20A region of the RGS data. First and second order extractions are plotted separately (m=-l in black, m=-2 in red) to 
facilitate line identification. The data from the two spectrometers have been averaged for each order extraction. The higher spectral resolution 
and resilience to source extent in is clearly seen in m=— 2 , where some line complexes blended in m =— 1 are resolved. 



Fig. 3. Illustration of the line 
profile used for estimating line 
fluxes. The plot compares the 
point source line spread func- 
tion for RGS1, the approxi- 
mate monochromatic line pro- 
file based on E0102’s angular 
distribution (derived from the 
Chandra image), and an heuris- 
tic wavelength broadening tem- 
plate function that is applied in 
addition to the angular distribu- 
tion. 


FWHM. The image was collapsed along the dispersion axis 
according to the roll of XMM-Newton during the observation. 
This 1-dimensional distribution was used as input to the re- 
sponse generator. Comparing the prediction to the brightest, 
isolated emission line (O VIII Lya), we see that the prediction 
is far narrower than the data, indicating a velocity distribution 
resolved by the RGS (see Fig. 3). Further accurate modeling 
requires coupling a line of sight velocity distribution to the an- 
gular distribution, but for the time being we are able to place 
a lower limit to the ejecta’s expansion, 1350 km s 1 f using the 
square velocity profile model of an expanding shell. This lower 
limit is consistent with the broadening seen in all of the lines. 
To proceed with spectral modeling, we use this distribution as 
a template for line modeling and flux estimation. 

We attempt to remove the effects of intervening absorption 
from the line intensities. The self consistent absorption column 
used is derived by fitting the long wavelength part of the spec- 


Table 1. Absorption corrected emission line ratios for single ion se- 
ries, for IE 0102-72.3. These ratios include resonant transitions only. 
The value used for Nh was 1.6 x 10 21 cm 2 , using solar abundances. 
Naturally, these ratios depend on the value of Nh assumed, particu- 
larly for C vi and O VII. Figure 4 plots these ratios along with pre- 
dicted ratios for CIE. Sample charge exchange (CE) measurements are 
included for O vm Lyman line ratios (Bliman et al., 1985), provided 
here for the discussion below. 


line ratio 

Carbon 

Oxygen 

Neon 

He a (nO/He /? 

NA 

6.3 ± .25 

7.9 ± .47 

(CIE) 


(9.2) 

(8.9) 

He /3/He 7 

NA 

2.6 ± .23 

2.2 ± .18 

(CIE) 


(3.2) 

(3.1) 

Ly a/Ly P 

5.2 ± .52 

6.8 ± .27 

7.4 ± .37 

(CIE) 

(12.5) 

(11-3) 

( 10 . 6 ) 

(CE, 0 +8 + H 2 ) 


(4.2) 


Ly p(Ly 7 

NA 

2.8 ±.17 

2.9 ± .29 

(CIE) 


(4.4) 

(4.1) 

(CE, 0 +8 + H 2 ) 


(3.1) 


Ly a/He a (w) 

NA 

1.35 ± .04 

0.84 ± .03 


trum (see the slope of the continuum in Fig. 1) with an ab- 
sorbed bremsstrahlung component hot enough to produce the 
Oxygen lines. The derived column (Nh = 1.6 x 10 21 cm 2 ) 
is twice the Galactic column toward the SMC and other esti- 
mates of absorption column toward E0102 (Blair et al., 1989; 
Hayashi et al., 1994; Blair et al., 2000). Absorption corrected 
line flux ratios are given in Figure 4 for multiple lines produced 
by the same ion. The line ratios inferred from the data are far 
from predictions for a collisionally excited, optically thin ther- 
mal plasma. As shown in the figure, the predicted intraseries 
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Fig. 4. Emitted line ratios. Plotted are the absorption corrected line 
ratios for adjacent emission lines for individual ions from Table 1, 
compared to line ratios expected from thermal plasma models. The 
left hand plot shows the ratios of Lya/Ly/? (C vi, O vm & Ne X) and 
Hec*(r)/He/? (Ovn & NelX); the right hand plot gives Ly/3/Lyi 
(OviII & Ne X) and He/?/He 7 (Ovn & NelX). The measured ratios 
are plotted as horizontal lines with error bars and the curves show the 
temperature dependence to these line ratios, assuming a CIE model. If 
contributions from recombination were turned off (the case of a purely 
ionizing plasma or strong NE1, for example), the curves would move 
up in value, particularly at the high T e end. 



Fig. 5. Derivation of the R (density) and G (temperature) line ratio 
diagnostics for O VII (left) and Ne IX (right). The data were fit using 
three line features, where we have used the heuristic wavelength dis- 
tribution function derived from O vm Ly a as a template to describe 
each (see Fig. 3). One parameter adjusted the ensemble line position 
(relative line positions were fixed), and the three line intensities were 
allowed to vary to fit the data. Individual contributions from the /, t 
and r lines are shown. The \l values are 1.4 and 0.9 for the O vn and 
Ne IX triplets, respectively. 


ratios a//? and / 3/7 for heliumlike and hydrogenlike ions are 
dependent on T et but asymptotically approach their minimum 
values in the limit of high temperature. Our measured ratios are 
significantly below these high temperature limits, for every ion 
species. It is curious to note that the measured ratios are similar 
to the expected line ratios from a recombining plasma. 

In spite of the angular extent of the source and the veloc- 
ity broadening of the lines, we resolve the Hea line complexes 
of O VII, Ne IX and Mg XI sufficiently to estimate contributions 


from the forbidden, intercombination and resonance transitions 
that provide robust diagnostics for the temperature and den- 
sity of the emitting medium (Gabriel and Jordan, 1969; Prad- 
han, 1982; Porquet and Dubau, 2000). Figure 5 shows the he- 
liumlike triplets of O VII and Ne IX, modelled using the known 
line wavelengths, and the line velocity profile adopted above 
(Fig. 3). While the density diagnostic R ratios obtained were 
reasonable and consistent with the low density limit at the 90% 
confidence level, the temperature diagnostic G values are low. 
They are measured at 0.55 ± .03 and 0.49 ± .05 for O VII and 
NelX, respectively. The electron temperatures corresponding 
to these ratios are 0.35 and 0.65 keV for O VII and NelX re- 
spectively (Pradhan, 1982), significantly above the maximum 
line producing temperatures (T m ) for each ion. These high tem- 
peratures are far from equilibrium temperatures where these 
ions are abundant, and suggest a purely ionizing plasma where 
the plasma is underionized and recombination rates are sup- 
pressed by the high T e . 


4. Discussion 

The Fe L lines that we do see in the SNR most likely trace the 
swept up ISM rather than the ejecta of the precursor. The Fe:0 
abundance ratio that we measure is about 1:900, compared to 
1:15 or so in the SMC (Russell and Dopita, 1992). With the ex- 
ception of the Fe, the composition of the SNR is dominated by 
products of He burning, the SN ejecta. The absence of emis- 
sion from N is suggestive that any nitrogen in the progenitor 
had been either blown off or burned into other products. 

The measured values for intraseries, single ion line flux 
ratios a//? and /?/7 are unexpected. They are not consistent 
with any T e capable of collisionally exciting the lines seen, nor 
any combination of such distributions. They are also far from 
values predicted from previous modeling efforts for adiabatic 
SNRs evolving through nonequilibrium ionization (Hamilton 
et al., 1983). 

We have been unable to identify any likely mechanisms 
by which an emitted electron impact excitation spectrum could 
be modified by selective absorption to produce the line ratios 
measured. Invoking an additional absorber fails: an intervening 
neutral absorber, used to explain the observed flux ratio pre- 
dicts columns different by Nh ~ 10 21 cm -2 between O VII 
and O VIII, which is clearly implausible. The possibility of the 
embedded gas working as an absorber is easily rejected: if sig- 
nificant, sharp increases in the Lya/Ly/3 ratios for Ovm and 
Nex, are expected due to the K edge energies of Ovn and 
Ne IX, which are not seen. 

The possibility is remote that resonant line scattering in the 
presence of a dusty or photoelectrically absorbing medium al- 
lows for the selective destruction of the large oscillator strength 
lines. In the approximation of total frequency redistribution, 
resonantly scattered photons are reemitted according to the doppler 
distribution of the gas, and the photon is likely to escape in r 0 
scatterings (where To is the optical depth at the line peak). The 
probability of being absorbed between scatters is r a fr # /ro, so 
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the fraction absorbed prior to escape is approximately r a bs * in- 
dependent of line oscillator strength. Even in the case of total 
coherence, where escape is likely after Tq scatters, the absorbed 
fraction is T a i 9 ro and we would expect the discrepancies in line 
ratios 6 * to scale with the ratios of oscillator strengths, which is 
not seen. 

As mentioned above, the line ratios, with the exception of 
the G ratios obtained for heliumlike ions, are reminiscent of 
a recombining plasma. Similar line ratios are also expected 
from charge exchange (Wise and Sarazin, 1989), where the ex- 
changed electrons are injected into specific high (n,/) levels 
according to energy conservation and details of the transient va- 
lence bonds formed between exchange partners. Spectra from 
such reactions have been observed extensively in tokamak and 
ion beam experiments, for example, in highly ionized Ar+ 18 
(Rice et al., 1986) and 0+ 8 (Bliman et al., 1985). In the lat- 
ter work, charge exchange reactions of 0 +8 with molecular H 2 
yielded values of Lya/Ly/3 ~ 4.2 and Ly/?/Ly 7 ~ 3.1 (40% 
relative error in each) and a preferential injection into n = 5. 
Calculations performed for interactions with atomic hydrogen 
yield values with comparable variances, but all well below the 
ratios predicted by electron impact excitation. 

In order for charge transfer to contribute significantly to the 
X-ray line emission from a SNR, low ionization species (neu- 
tral hydrogen or helium, or singly ionized helium) must come 
into close proximity with more heavily ionized metals (i.e., 
fully stripped carbon, oxygen, and neon, for this case). Given 
that the charge exchange and electron impact excitation rate co- 
efficients are of order 7 ce ~ 10“ 6 andTEiE ~ 10" 11 cm 3 s _1 , 
respectively, charge exchange can be important if the neutral 
fraction for hydrogen is greater than ~ 10~ 5 or so (Wise and 
Sarazin, 1989). Unfortunately, this is well above the neutral 
fraction expected for the hot ionized plasma where the X-ray 
emission is thought to originate. 

However, there is abundant evidence that the X-ray emit- 
ting media in SNRs are far from homogenious. Cool clouds 
can coexist in rough pressure balance with the hotter inter- 
cloud medium. One therefore expects a steady flux of highly 
ionized material from the hot medium at the “faces” of each of 
the clouds. These ions will charge exchange with neutral hy- 
drogen in the clouds, giving rise to X-ray emission lines. Since 
the charge exchange cross-section for the ions is considerably 
larger than the photoelectric absorption cross-section for the 
emitted line photons, the emitted line radiation should emerge 
and contribute to the observed X-ray spectrum of the remnant. 
This process should be important if the mean intercloud sepa- 
ration is smaller than the characteristic distance d an ion will 
move through the hot gas before it is collisionally excited: 

d- l.Opc (7£/£/10 -u ) 1 W 

The Chandra image of E0102 (Gaetz et al., 2000) exhibits 
significant structure on spatial scales comparable to or less than 
this value. It therefore seems likely that charge exchange at 

6 The effect on a//? ratios should be about twice the effect on /?/7 

ratios. 


cloud interfaces plays a significant, if not dominant role in pro- 
ducing the anomalous line ratios that we observe. 
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